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The zero-energy bound states at the edges or vortex cores of chiral p-wave superconductors should 
behave like Majorana fermions. We introduce a model Hamiltonian that describes the tunnelling 
process when electrons are injected into such states. Using a non-equilibrium Green function formal- 
ism, we find exact analytic expressions for the tunnelling current and noise and identify experimental 
signatures of the Majorana nature of the bound states to be found in the shot noise. We discuss the 
results in the context of different candidate materials that support triplet superconductivity. Exper- 
imental verification of the Majorana character of midgap states would have important implications 
for the prospects of topological quantum computation. 



When the individual constituents of a many-body sys- 
tem interact non-trivially with each other, they can give 
rise to low-energy states in which the elementary ex- 
citations are very different from the original building 
blocks. Examples among electronic condensed-matter 
systems include the spinous and holons in Luttinger liq- 
uids (realized, e.g., in single-wall carbon nanotubes) or 
the Laughlin quasiparticles of fractional quantum Hall 
systems. In the context of superconductors, Cooper pairs 
can be seen as a relatively simple example of such ex- 
citations, but more exotic states are also possible. In 
the present work we are interested in the case of p-wave 
chiral superconductors (i.e. with an order parameter of 
the type ± ipy), examples of which would be stron- 
tium ruthenate (Sr2Ru04) [l|, Q and, possibly, a num- 
ber of organic superconductors like the Bechgaard salts 
({TMTSFjsX; X=PF6,C104,...) 0, and even heavy 
fermions {e.g. UPta) 

Superconductors with p-wave orbital symmetry have 
spin-triplet pairing and the order parameter is a tensor 
in spin space rather than a scalar. This introduces ex- 
tra freedom and allows for different types of supercon- 
ducting phases, first studied and observed in superfluid 
■^He. In the so called A-phase, Cooper pairs are in a 
state dubbed 'equal (or parallel) spin pairing' (ESP); all 
the examples given above are candidate systems for ESP. 
Within weak-coupling BCS theory the up- and down-spin 
sectors are then independent from each other and the re- 
spective Bogoliubov-de Gennes (BdG) equations are de- 
coupled. Another aspect of the A-phase of p-wave su- 
perconductors is that it can support vortex-core bound 
states with a spectrum given by — — + with 
n (z 1j, n — and loq a frequency that depends on the 
details of the vortex profile ■&]. In particular, for n = 0, 
one notices that the vortices support 'zero modes'. (This 
should be contrasted with the case of s-wave supercon- 
ductors for which the vortex bound-state spectrum has 
again the same form but this time n — 1/2, the other of 
the only two possibilities consistent with the E ^ —E 
symmetry of the BdG equations.) A more detailed con- 
sideration of these midgap bound states in the case of 
ESP reveals that they have self-adjoint wavefunctions. 



naturally described as Major aria fermion modes, and can 
also be found as edge states i6!|. Such Majorana states 
constitute one more example of an exotic low-energy col- 
lective excitation that is very different from the original 
electrons that condensed into the superconducting state. 

It would be already extremely interesting to be able to 
experimentally identify these strange Majorana bound 
states, since that would constitute a stringent test of our 
current picture of EPS p-wave superconductivity, but 
the motivations run further. The availability of Majo- 
rana fermions can be exploited in the context of quan- 
tum computation, a completely new and revolutionary 
approach to computing that would mix aspects of the 
digital and the analog computing paradigms by exploit- 
ing the basic laws of quantum mechanics. Majorana op- 
erators (call them 7]i) can be taken in pairs to define stan- 
dard fermionic operators [say, = [rj^ + irji) /"v/S]; each 
of these generates a two-dimensional Hilbert space that 
can be used to define a quantum-bit (qubit). Because 
the two Majoranas can be spatially far apart {e.g. in 
two different vortices) and are very different from the 
usual fermionic quasiparticles around them, the so de- 
fined qubit would be relatively immune to decoherence 
,7;], which would sidestep one of the crucial problems 
faced by the development of quantum computing hard- 
ware. Moreover, it turns out that the usual global gauge 
symmetry of the fermi fields is reduced to a discrete Z2 
symmetry for the Majoranas at the core of a vortex (and 
they can be shown to change sign when a third vortex 
moves about encircling them 

UM)- Chan ging the sign 
of a single Majorana of the pair that defines a qubit op- 
erates the change ^ c, or, in other words, it acts as 
a qubit-flip (or q-NOT) gate, and the Z2 symmetry be- 
ing discrete leaves no room for errors. This shows that 
braiding vortices would perform quantum-logical oper- 
ations on the information stored in them, an approach 
known as topological quantum computation (for a recent 
review see Ref. 9). 

Recently, the overlap matrix element between a local- 
ized electron and a Majorana bound state was computed 
for a model of a superconducting wire (introduced in the 
context of quantum computation |l Cl| ) and found to be 
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non-zero This indicates that tunneUing transport 

into Majorana modes is possible, which opens interesting 
possibilities since tunnelling has proved repeated times 
to be an invaluable tool in the study of superconducting 
states. The study of tunnelling noise might also be use- 
ful, since shot noise is another probe able to distinguish 
normal versus superconducting states and ballistic ver- 
sus diffusive transport l3, 13] ; it can also be sensitive to 
the charge and statistics of the carriers and was used, for 
instance, to confirm the presence of Laughlin quasiparti- 
cles in fractional quantum Hall devices |14| . Noise probes 
can be local, in order to study localized states [15[ . The 
purpose of this article is to model the tunnelling pro- 
cesses into Majorana bound states and to determine the 
current and noise characteristics in order to identify sig- 
natures that would allow the experimental identification 
and study of such states. A generalized geometry of the 
experiment we consider is shown in Fig. [TJ We shall find 
that the Fano factor (or shot noise to current ratio) for 
such tunnelling processes has unit matrix structure and 
is given by 



Sap [UJ = 0) 
y/T-^oo e {la + I 13) 



lim 
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where a, /3 = {L, i?} = ±1 label the lead where the cur- 
rents {la, 13) are measured and Safi is the noise spectrum 
defined below. This is different from the result for a regu- 
lar fermionic bound state for which FajS — 1/2 has a 'fiat' 
matrix structure and is half as big -the full expressions 
for the noise are given in Eqs. ^ and ([3])-. We shall 
argue that measuring the Fano factor would provide a 
clear signature of the Majorana nature of a bound state. 

The current-voltage characteristics for tunnelling into 
low-dimensional chiral p-wave superconductors was com- 
puted before for the case of planar junctions using a 
BTK formalism (including the bound states in a den- 
sity of states approximation) [l^ . or for point contacts 
using a microscopic tunnelling Hamiltonian and non- 
equilibrium Green functions but in the absence of bound 
states {i.e. far from vortices or edges) [l3|- Here we 
concentrate on the tunneling into bound states or edge 
modes for voltages smaller than the superconducting gap 
and carry out full microscopic calculations for the current 
as well as the noise. Our starting point is the following 
tunnelling Hamiltonian: 

H = Ho + V2tL im + Vl) + V2tR - t/.^) rjr 

here r]r,i are two different Majorana operators (located, 
for instance, in two different vortices; see Fig. [T|) and 
V'i = / V'lfeC^fc/27r are the fermions at the position of the 
point contact in each of two leads. We consider only the 
relevant spin projection of the ESP state and effectively 
work with spinless fermions. The Hq term in the Hamil- 
tonian has two parts. 
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FIG. 1: Schematic representation of our model for a setup in 
which two contacts {e.g. two STM tips) tunnel electrons into 
the Majorana bound states at the core of vortices. The two 
Majorana states might be linked with a tunnelling amplitude 
e. The setup is generic and any or both of the vortices could 
be replaced by edge modes or obviated altogether. The model 
would also apply to the case of multiband vortices with two 
different Majoranas one on each band, in such a case e corre- 
sponds to the amplitude of band mixing that can take place 
at the core of the vortex. 



where is defined as above and there is a term i?o,a 
for each lead. i/o,c serves to model the case when the 
overlap matrix element between the two Majorana bound 
states is non-zero (cf. Ref. Illh. Notice that the terms in 
H involving 77 's are bound to have the form they have 
due to hermiticity requirements (cf . Ref. @) . The overlap 
amplitudes ta have to be real and we can take them to 
be positive for the sake of concreteness. The choice of 
relative signs in the tunnelling terms is arbitrary and 
amounts to a choice of global gauges for the leads. 
For generality, we shall rewrite the Hamiltonian as 



[ta {c^'lpa +1pic) + aSa {iplc^ + Clpa)] 



a=L,R 



e irii'qr 



+e c'c 



The Majorana tunnelling case is recovered by setting 
5a ~ ta, but on the other hand we can set Jq. = and we 
have the standard resonant level model that would de- 
scribe a regular {i.e. non-Majorana) edge or bound state 
[ist . We can thus discuss the two cases using a common 
familiar language and compare them more easily. 

We start by computing the current, which is given by 

I{t)^^ {dt {Nn - Nl)) = ^.{[H,Nl- Nn]) 
2 2,1 

a=L,R 

where is the normal Keldysh component of the equal- 
time Green function, iG^^ = (V'ctC''')j^^j, and is its 
anomalous counterpart, iF^^ — {^aC)^^\ (and similarly, 
mutatis mutandis, for other components). In order to 
compute these Green functions we follow the 'local action 
approach' of Ref. [iSj, with the main difference that here 
the calculation is done fully analytically. The case of 5a — 
is relatively straight forward and was discussed before 
in the literature (see pil for a review); we therefore give 
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here only the results for the Green functions when 6a = 
ta- In order to compute them we work with positive 
frequencies and adopt the spinor basis given by 

= vt(^) ^^(^) ^t^(^) ^(^) ^t(^))^ 

(where the bars stand for minus signs). Inverting the 
local action -which is equivalent to solving the full non- 
equilibrium Dyson equations- and restricting ourselves 
to the symmetric case (II = tj^ ~ t), we obtain the 
following Green functions: (i) the localized-states Green 
fimctions, 

G-4f.De(c.-e)(th^^+th^ 

e'4{-7^cH (th^ + th^ 

where we used T = + Tj^ with = ^ {ti + 5l) 
(later we will need also T = ^ OiSata) to define the 
functions {uj) = -^J^ and {u) = ;;j2r^-f)7- The 
sums on 7 = {L, R\ are implicit and W — Avy (wf is 
the fermi velocity in the leads). We introduce now the 
notation Lp — —La — 1 and write, (ii), the inter/intra 
lead Green functions, 



We now use the third set of Green functions and re- 
place in the formula for the current. Let us first quote the 
result for a resonant level or double barrier (cf . Ref . Il3l ) , 
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a=R,L 
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Notice the current becomes zero if either V ^ or vanish; 
in fact, we always find II — Ir (with la = ea (dtNa))- 
The situation is different when the tunnelling is into a 
Majorana state; in such a case, the two currents are in- 
dependent (parametrically related if e 7^ 0) and given by 
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with (valid also if =/= tji) 



Dca (t^) 



X [Dc {LJ-e)+ pjF* {uj) + a-fFc (lj) + a(3Dc (w -I- e)] 
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th ■ 
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X [aDc {oJ-e) + aPjF^ {lu) + jF^ {uj) + (30^ {oj + e)] 
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th 
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Finally we use the convention that upper/lower indices 
correspond to upper/lower signs (>w= means there is no 
complex conjugation) and write, (iii), the 'tunnelling' 
Green functions. 
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If Tl = F/j and e = 0, then / = {II + Ir) /2 coincides 
with the result for a resonant level. Even though the 
results are mathematically different (notice that for the 
case of a wire, e 7^ always regardless of its size (llj), 
the differences might be hard to detect experimentally; 
but we shall now see that the noise provides a more ro- 
bust signature of the Majorana nature of the tunnelling 
intermediate states. 

The noise is a measure of the deviations of the current 
from its average value [Ala (t) — la (t) ~ {la (t))] and it 
is standard to define it as the following correlator |13[: 

Sap {t (A/o (<) Alp it') + Alp (f) Ala (t)) . 

Its Fourier transform is known as the noise power spec- 
trum and for i5q = is given by 
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For the general case, the expression is similar but longer; 
it involves also anomalous Green functions and comprises 
thirty two terms with a mix of both correlation 2J| and 
convolution products. 

We concentrate on the zero- frequency noise component 
and first rederive the result for a resonant level 13, 20l : 
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where eV = fiL — fJ.R- In the fully-symmetric case {fiL 
~^J'R = M ^-iid = tfi = t) it takes the form 




Sep (0) 



Lu — e 



r2 



(2) 

and the Fano factor becomes Fap = 1 — 2TRrL/^^ 
1/2 as mentioned earlier. On the other hand, when the 
tunnelling takes place into a Majorana bound state we 
find the following result (for the fully-symmetric case): 



Sa0 (W = 0) 



= cth 



eV_ 



25al3- 



i-^s (e - ijf + r2 



Notice that the diagonal and off-diagonal matrix com- 
ponents of Sai3 are different now. In particular, we re- 
mark that lime_,o Saa = 0. Taken together with the 
result given above for the current, this indicates that in 
the e — > limit the right and left tunnelling processes 
are completely independent even at the level of current 
fluctuations. It is therefore instructive and important to 
study this case more in detail, because of its greater sim- 
plicity and its relevance to single-tip setups. We relax the 
condition on the chemical potentials {i.e. consider fi^, fiR 
arbitrary) and find that the noise can be written as 



Sep (0) 
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Given the right-left independence, we expect the expres- 
sion to be valid also when II tn. In particular, this im- 
plies that the Fano factor, Eq. ([1]), is not sensitive to the 
contacts asymmetry, unlike what happens for (5_r,l = 0. 

Presently, the likely best place to find single-Majorana 
bound states is at the edges of Bechgaard salts. In the 
presence of magnetic fields much larger than the param- 
agnetic limit, but still smaller than Hc2, the supercon- 
ducting state ought to be similar to the Al-phase of '^He 
in which all the spins are aligned with the magnetic field 
and effectively there is only one spin species. The situa- 
tion in the case of strontium ruthenate is more complex 
and the existence of single-Majorana bound states seems 
more elusive. Let us start by pointing out that the vor- 
tices that support isolated Majorana zero modes are no 
less bizarre themselves: they carry only half of the super- 
conducting flux quantum and the vorticity lies entirely in 
one of the spin sectors while the other one does not show a 
winding phase at all 21| . Recent scanning tunnelling mi- 
croscope (STM) studies of Sr2Ru04 have found a square 
lattice of vortices with a full flux quantum each , even 
though the magnetic field was presumably large enough 
to take the system into an ESP state [l|, lof. However, 



the experiments did find a strong zero bias conductance 
peak that remains unexplained. One possibility is that 
the observed vortices are built out of two half-vortices, 
one for each spin projection. In that case, one would 
expect, respectively, two Majorana bound states, and e 
in our model would be related to the amplitude of spin 
mixing at the vortex core. If this is the case, the vor- 
tices would not have the same braiding properties as the 
half-vortices, but the Fano factor would nevertheless be 
sensitive to the Majorana nature of the midgap states. 
Whether these 'double half-vortices' would still be useful 
for topological quantum computation would require fur- 
ther investigation, but identifying them experimentally 
would be extremely interesting in any case. 
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